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O ; ABSTRACT 
, 

We present three-dimensional MHD simulations of buoyant magnetic flux tubes that 
. rise through a stratified model convection zone in the presence of solar rotation. The 

l/^ , equations of MHD are solved in the anelastic approximation, and the results are used 

I to determine the effects of solar rotation on the dynamic evolution an 0-loop. We find 

I that the Coriolis force significantly suppresses the degree of fragmentation at the apex 

' of the loop during its ascent toward the photosphere. If the initial axial field strength of 

• the tube is reduced, then, in the absence of forces due to convective motions, the degree 

of apex fragmentation is also reduced. We show that the Coriolis force slows the rise 



of the tube, and induces a retrograde flow in both the magnetized and unmagnetized 



c/2 , plasma of an emerging active region. Observationally, we predict that this flow will 

^ ' appear to originate at the leading polarity, and will terminate at the trailing polarity. 

^ ' Subject headings: methods: numerical — MHD, Sun: interior. Sun: magnetic flelds 



1. Introduction 



Active regions represent the largest observable concentrations of magnetic flux on the Sun. 
The magnetic fleld within active regions exhibits a bipolar structure, which suggests that they are 
the tops of large O-shaped loops which have risen through the convection zone and have emerged 
through the photosphere. On average, active region bipoles are oriented nearly parallel to the 
East- West direction (Hale's Polarity Law). This suggests that the subsurface field geometry is 
approximately toroidal. Since active regions appear to be the tops of large, rising Q-loops, and since 
Hale's law persists for years in a given solar cycle, we can infer that the toroidal flux must exist well 
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below the solar surface, and that the magnetic flux must reside in a region where it is relatively 
free from disruption by buoyant instabilities or convective motions (Spruit & van Ballegooijen 
1982; van Ballegooijen 1982; Ferriz-Mas k, Schiisslcr 1993, 1995). However, if the magnetic flux is 
embedded in the radiative zone (well below the convective zone) it would be so stable that it could 
not emerge through the surface on the time scale of a solar cycle (Parker 1979; van Ballegooijen 
1982). This implies that the toroidal fleld likely resides in the "convective overshoot region" — a 
thin, slightly convectively stable layer located at the interface of the radiative and convective zones. 
This layer also appears to coincide with the "tachocline" (Kosovichev 1996; Corbard et al. 1999), 
where the solid body rotational behavior of the radiative interior transitions into the differential 
rotation behavior observed in the convection zone. Thus, the toroidal fleld seems to be generated 
and stored in the overshoot layer. (Oilman, Morrow, & DeLuca 1989; DeLuca & Oilman 1991; 
Parker 1993; MacGregor & Charbonneau 1997; Durney 1997; Dikpati & Charbonneau 1999). 

We now have a basic picture of how active regions might form: A portion of the toroidal flux 
layer eventually succumbs to one of several possible instabilities (eg. Caligari et al. 1995; Fan & 
Fisher 1996; Fan 2000; Wissink et al. 2000b), and forms a magnetically buoyant tube that begins to 
rise through the convection zone toward the photosphere. The top of the emerging flux loop is then 
observed at the photosphere as a bipolar magnetic region. The simplest way to model this picture 
is with the "thin flux tube" approximation (see Defouw 1976; Roberts &: Webb 1978; Spruit 1981). 
In this simple model, it is assumed that magnetic flux tubes move through a field-free plasma, that 
the cross-section of the tube is much smaller than all other relevant length scales of the problem, 
and that pressure balance is always maintained across the tube. Given these assumptions, one 
can derive an equation of motion for a one-dimensional tube moving through a three dimensional 
medium. This model has had great success in explaining many of the observed properties of active 
regions, including (but not limited to) the variation and scatter of active region tilt as a function 
of solar latitude (D'Silva & Howard 1993; D'Silva & Choudhuri 1993; Fan et al. 1993; Galigari et 
al. 1995; Fisher et al. 1995; Longcope &; Fisher 1996; Fan & Fisher 1996), the asymmetries in the 
orientation of active regions just after emergence (van Driel-Oesztelyi & Petrovay 1990; Moreno- 
Insertis et al. 1994), and the latitudinal variation and fluctuations of fleld line twist from active 
region to active region (Longcope et al. 1998). 

However, from the start, the thin flux tube approximation assumes that the magnetic field will 
retain its "tube-like" identity; that is, it will not fragment or disperse as it evolves. This assumption 
is called into question by recent two-dimensional MHD simulations which show that unless flux tubes 
possess a critical amount of initial field line twist along the tube — and observations of emerging 
active regions suggest that many do not (Pevtsov et al. 1995; Longcope et al. 1998, 1999) — the 
tubes will fragment (break apart) before they are able to emerge through the surface (Linton et al. 
1996; Moreno- Insertis &: Emonet 1996; Emonet &: Moreno-Insertis 1998; Fan et al. 1998). Recently, 
Abbett et al. (2000) used three-dimensional MHD simulations to demonstrate that the degree of 
cohesion at the apex of an emerging fi-loop depends strongly on the three-dimensional geometry 
of the loop. If the curvature of the J7-loop at its apex is relatively large (ie. if the loop is tall and 
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narrow), then the apex will emerge through the photosphere more cohesively than if the curvature 
is small (ie. if the J7-loop is short and wide). Dorch & Nordlund (1998) and Dorch et al. (1999) 
also performed detailed three-dimensional MHD simulations of buoyant magnetic flux tubes. All 
of these 3-D studies still require a significant amount of initial field line twist in order to suppress 
the eventual fragmentation of the tube near the loop apex. Is there a way out this quandary which 
Abbett et al. (2000) refer to as "Longcope's paradox" ? 

In this paper, we investigate how solar rotation affects the dynamic evolution of buoyant 
magnetic flux tubes. We find that the Coriolis force plays a dynamically important role during 
the tube's ascent, and acts to maintain the tube's cohesion even in the absence of initial field line 
twist, possibly resolving Longcope's paradox. In Section 2 we briefiy describe our computational 
methodology, discuss how we include the Coriolis force in our computational model, and define the 
relevant physical parameters of our problem. In Section 3 we analyze the results of our simulations, 
and develop a very simple model that effectively characterizes the simulation results pertaining to 
flux tube cohesion. In this section we also comment on possible observational consequences of the 
simulations. Finally, we present our summary and conclusions in Section 4. 
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2. 



Method 



We use the methods described in Abbett et al. (2000) and Fan et al. (1999) to solve the three- 
dimensional equations of MHD in the anelastic approximation (Ogura k, Philhps 1962; Gough 
1969; Lantz k, Fan 1999). The anelastic equations result from an expansion of the equations of 
compressible MHD about a zeroth order reference state (which we take to be a ficld-frcc, adiabati- 
cally stratified polytrope of index m = 1.5). In this approximation, the higher order time derivative 
of the perturbed density is neglected, effectively filtering out fast-moving acoustic waves. This is 
computationally advantageous, as it allows for much larger timesteps than would be possible in a 
fully compressible calculation. This approximation is valid in sub-surface regions where the acous- 
tic Mach number is small (M = v/c^ <^ 1), and where the local sound speed (cg) greatly exceeds 
the Alfven speed {va) of the plasma (Glatzmeier 1984). A detailed theoretical discussion of the 
anelastic formulation can be found in Lantz k, Fan 1999 (and references therein). 

The anelastic equations in a non-rotating reference frame are given by equations [l]-[7] in 
Abbett et al. (2000) (hereafter AFF). To consider the effects of solar rotation, we adjust the anelastic 
momentum equation to include non-inertial terms (we note, here, that due to a typesetting error in 
the final version of AFF, a cross-product symbol was omitted in the Lorentz force term of equation 
[2]). To consider the effects of solar rotation, we adjust the anelastic momentum equation to include 
non-inertial terms. In this study, we neglect centrifugal effects and concentrate only on the effects 
of the Coriolis force due to solid-body rotation. In this approximation, the anelastic momentum 
equation (in the rotating reference frame) becomes: 



Here, po refers to the density stratification of the zeroth order reference state, while pi, pi, v, and B 

represent the density, pressure, velocity and magnetic field perturbations. The viscous stress tensor 
is given by = fi{dvi/dxj + dvj/dxi — 2/3(V • v)(5jj), where fi is the coefficient of viscosity and 5ij 
is the Kronccker delta function, ft denotes the angular velocity about the Sun's axis of rotation. 
Both the acceleration due to gravity, g = — g z, and the coefficient of viscosity are assumed constant 
in our calculations. 

We solve the non-dimensional form of the anelastic equations in the modified, local /-plane 
approximation of Brummell et al. (1996). This approach is similar to the standard /-plane approxi- 
mation of geophysical fluid dynamics (Pedlosky 1987), except that the transverse components of the 
Coriolis force within the rectangular domain are retained. Since is assumed constant throughout 
the Cartesian box (positioned at a given latitude), we can apply simple, periodic boundary con- 
ditions in the horizontal directions, and stress-free, non-penetrating conditions at the upper and 
lower boundaries. Further simplification is possible by noting that the momentum density (in the 
anelastic approximation) and magnetic field vectors are both divergence-free. This allows us to 
define scalar potentials such that B = V x V x Bz -I- V x J'z and 

PqV = V X V X Wz -I- V X Zz . Along with the other dependent variables of the system, these 




) 



Vpi Pig - 2po (fl X v) -h — (V X B) X B + V • n . 



(1) 
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potentials are decomposed spectrally in the horizontal directions: 

Q = i;Qmn(^)e''^^(^-^ + ^"J'). (2) 
mn 

Here, Q and Qmn refer to a given dependent variable or scalar potential and its associated Fourier 
coefficient, /m and are defined as m/>C-r and n/£y respectively, where and Cy denote the size 
of the domain in the x and y directions. The integers m and n range from — N-r/2 + 1 to Nj;/2 
and —Ny/2 + 1 to Nx and Ny refer to the number of zones in the x and y directions. The 

Fourier coefficients are then discretized with respect to the vertical direction (defined to be parallel 
to g), and vertical derivatives are approximated by fourth-order, centered finite differences. Using 
the semi-implicit method of operator splitting, we apply the second-order Adams-Bashforth scheme 
to the advection terms, and the second-order Crank-Nicholson scheme to the diffusion terms (see 
Press et al. 1986 for a description of these methods). 

A detailed, in-depth description of the techniques we employ to obtain a solution to the anelas- 
tic equations in a non-rotating reference frame can be found in Appendix A of Fan et al. (1999) 
(hereafter FZLF). Only minor modifications are necessary to extend that formalism to apply in 
the case of a rotating reference frame in the modified /-plane approximation. The equation for 
the time evolution of Z^an, the Fourier coefficients of the spectral decomposition of Z, is obtained 
by evaluating the i-component of the curl of the momentum equation. Thus, the addition of the 
Coriolis term to equation (1) results in an additional non-inertial term to equation (A16) of FZLF: 

{ ^ ) = + ATTiW^n iUn. + gn^y) . (3) 

V / cor 

Here, the Jlj's represent the Cartesian components of the angular velocity vector, and Wmn denotes 
the Fourier coefficients of the scalar potential W. Equation (A18) of FZLF results from taking 
the X and y components of the curl of the momentum equation as it appears in equation (A4). 
Thus, we must add another term to equation (A18) to reflect the addition of the Coriolis force to 
equation (1): 

\ ot J oz y po J poho Po 

In the above equation, = ~^(27r)~^(/^ -|- g^)~^ {fm^^ymn — 5n^^xmn) is the scalar variable that 
corresponds to Omn in equation (A18) of FZLF (we change notation so as to avoid confusion with 
the components of the angular velocity vector); ujxmn and tOymn are the Fourier coefficients of the 
X and y components of the vorticity; and /iq = —po{dpo/dz)~^ is the density scale height. The last 
necessary modification is to equations (A23) and (A24) of FZLF where we add the contribution from 
the horizontal average of the Coriolis term in the momentum equation: {d{poVx) / dt)cor = "^Po^zVy 
and {d{poVy) / dt)cor = —"^Po^z^x- Here, Vx and Vy denote the horizontal averages of the x and y 
components of the velocity (note that Vz is zero). 
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We now turn our attention to the non-dimensional formulation of the anelastic equations as 
implemented in the code. As was described in AFF, many of the dimension-less physical variables 
and fundamental physical parameters are expressed in terms of the magnetic field strength along the 
axis of the initial flux tube, Bq. With the addition of the Coriolis force to the momentum equation, 
we introduce another fundamental dimension-less quantity, Rb = VA/{2^lHr) = 1/(20'), which 
we refer to as the "magnetic Rossby number" (a quantity also discussed by Schiissler & Solanki 
1992). Here, fi' is defined as the non-dimensional rotation rate given in terms of the ratio of the 
characteristic velocity scale to the characteristic length scale of our problem; namely, the pressure 
scale height at the base of the computational domain (Hr), and the Alfven speed along the axis of 
the initial tube {va = Bq/ where the characteristic density, p^j is taken to be the density 
at the lower boundary). The magnetic Rossby number is a measure of the relative importance 
of the Coriolis force on the dynamic evolution of our buoyant magnetic flux tube. Rotational 
effects become important when Ri, < 1, and are less significant when i?^ > 1. We therefore wish 
to investigate three representative cases: i?^ = 1/2, R^ = 1, and Ri, = 2. Of course, the solar 
rotation rate is not a free parameter — estimates of in different parts of the convection zone 
can be obtained from helioseismologic measurements. Thus, varying the magnetic Rossby number 
to investigate the effect of solar rotation on a buoyant magnetic flux tube amounts to varying the 
initial axial field strength of the tube. 

The specification of Bq impacts the scaling of two of the non-dimensional diffusive parameters: 
the Reynolds number [Re = [p][^]HAt~''^), and the magnetic Reynolds number [Rm = [^]H'7~^)- 
Again, the characteristic length, velocity, and density scales {[I], [v\, and [p]) arc given in terms 
of Hr, VA, and respectively. For this study, we hold the coefficients of viscosity and magnetic 
diffusion [p and r/ respectively) constant between runs (where possible), and vary Bq such that we 
obtain our three representative values of the magnetic Rossby number. In this manner, we focus 
our attention on the effects of solar rotation on the dynamic evolution of a rising fi-loop. We note 
however, that as Bq becomes large, it is increasingly difficult to hold the coefficient of viscosity 
constant for a given grid resolution. This is because as the Reynolds number increases, features 
develop at the forefront of the rising tube and in the trailing eddies that are on the order of a 
few grid zones in size. If this occurs, the code becomes unstable, and the grid resolution must be 
increased to alleviate the problem. Table 1 lists the values of the relevant physical parameters for 
each of seven simulations. For the cases where Bq ~ lO^G, values of p and r/ were increased to 
avoid instabilities. We otherwise would have to double the grid resolution in the tube cross-section 
to take full advantage of the FFT algorithms employed in the code, and such a large increase is 
unnecessary for the purposes of this study. 

Our labeling convention (sec column one of Table 1) consists of four letters which signify both 
the strength of Bq, and the latitude of the computational box. For example, LFHL refers to the 
run which has an initial magnetic flux tube with a relatively Law axial Field strength positioned 
at a i/igh Latitude, while HFLL refers to the run with an initial tube of relatively ifigh axial 
Field strength positioned at Low Latitude — and so on. Run SSO refers to the simulation where 
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solar rotation was neglected (this notation was chosen in order to be consistent with the naming 
convention of AFF). In each case, we begin with a static, horizontal, cylindrical magnetic flux tube 
positioned near the bottom of a simulation box that spans 5.147 pressure scale heights vertically. 

2/2 

The initial field configuration is given by B = Bx{r)x, where Bx{r) = BQe~^ and r is the radial 
distance to the central axis in the tube's cross-section. The size of the initial tube, a, is taken to 
be O.lHr in each of the simulations. We choose not to include any initial field line twist about the 
central axis of the tube, since we wish to determine what role the Coriolis force alone plays in the 
overall fragmentation of the tube as it rises toward the photosphere. As was done in AFF, we apply 
an artificial entropy perturbation to the midsection of the tube at time t = 0. This causes the tube 
to rise as an Q-loop without having to take the computationally expensive step of evolving each 
case self-consistently from a state of initial force balance (as was done in Caligari et al. 1995; Fan 
& Fisher 1996; Fan 2000). 



3. Results 

Each simulation begins with an untwisted cylindrical magnetic flux tube embedded near the 
base of our model convection zone. An entropy perturbation is applied to the flux tube such 
that in the absence of rotation, the tube will rise as an $7-loop, and the apex will approach the 
photospheric boundary in a highly "fragmented" state. That is, a majority of the magnetic flux 
will be concentrated away from the tube's central axis along vortex pairs that are shed during the 
tube's ascent. Following AFF, we consider a section of the O-loop to be fragmented if the ratio 
of the magnetic field weighted second moments of position along the Frenet binormal and normal 
directions of the tube exceed 1.5 (see Section 3 of AFF for details on how tube fragmentation is 
determined, and see Figures 1 and 3 in AFF for visual examples of a fragmented and non-fragmented 
Q-loop). We expect that solar rotation will play a dynamically important role during the tube's 
ascent, since the loop's rise time is on the order of a rotation period for reasonable initial field 
strengths (3 x lO^G < Bq < 1 x lO^G, see Fisher et al. 2000). 

We find that the impact of the Coriolis force on a rising flux tube is dramatic, as demonstrated 
in Figure 1. In this figure, we compare the magnetic field strength in the apex cross-section of two 
O-loops: one which rises at 15° latitude in the presence of solar rotation (J7 ~ 13.9° /day) with an 
initial axial field strength of Bq = 2.67 x 10^ G (corresponding to Rb = 0.5), and an identical loop 
which rises in the absence of rotation. Two effects are immediately apparent. First, the il-loop that 
rises through the rotating atmosphere has yet to attain the height of the loop that rises through the 
non-rotating atmosphere; and second, the loop subjected to the effects of the Coriolis force retains 
its cohesion, and does not fragment. Figure 2 is a volume rendering of a rising magnetic flux tube 
in a run where Rj, = 2 and Bq = 1.0 x 10^ G. The coordinate axes are chosen such that y points 
northward, and z points radially outward. Again, we see that the tube retains its cohesion during 
its ascent. Figure 2 also demonstrates that in the presence of rotation, the fi-loop is somewhat 
asymmetric about the apex. That is, the trailing portion of the loop forms a steeper angle with 
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respect to the x direction than the leading side of the loop. 

A qualitative understanding of these effects can be obtained by examining Figure 3, which 
shows |B|, the transverse components of the Coriolis force, and the flow field in the apex cross- 
section during the initial stages of run LFLL (an earlier snapshot of the same run that is displayed 
in the left frame of Figure 1). In this plot (and all others, unless otherwise stated) the physical 
quantities are displayed in naturalized units: the magnetic field strength is given in terms of the 
strength of the field along the axis of tube at t = {Bq), the velocities are given in terms of the 
Alfven speed along the axis of the initial tube {va), and the unit of time is the ratio of the pressure 
scale height at the base of the computational domain to the characteristic velocity, Hr/vA- For run 
LFLL, shown in Figure 3, Bq = 2.67 x 10^ G, va = L68 x 10"^ cm s-^ (for pr = 0.2 g cm-^), and 
t = 1.5 corresponds to 6.18 days. The axes in the figure are given in terms of computational zones; 
in this case, each individual zone is O.OlTifr, or approximately 1 Mm. 

In the right-hand frame of Figure 3 we sec that a strong, retrograde axial flow has developed 
as a result of the Coriolis force acting on rising material within the buoyant loop. This is in sharp 
contrast to the simulations in AFF, where solar rotation is neglected and the initial flux tube has 
little or no field line twist. In these cases, the strong axial flow is absent, and counter-rotating vortex 
pairs can quickly form on either side of the tube's central axis (see also Wissink et al. 2000a; Dorch 
&; Nordlund 1998). As was shown in AFF and Longcope et al. (1996), the non-vertical component 
of the hydrodynamic "lift" force per unit volume generated by flows in the vortex tubes acts to push 
them apart, ultimately resulting in fragmentation near the apex of the il-loop. However, in the 
presence of rotation, we find that the circulation (and thus the degree of fragmentation of the loop) 
is substantially reduced. Roughly speaking, this follows from the fact that the Coriolis force can 
do no work (it always acts in a direction perpendicular to the flow); it simply redistributes kinetic 
energy between transverse and axial flows. The initial increase in the axial flow at the apex of the 
loop gives rise to a secondary effect: a transverse component of the Coriolis force normal to the 
angular velocity vector. At low latitudes, this component acts primarily in the —z direction, and 
opposes the buoyant rise of the loop; however, it also forces the top of the loop to drift northward 
over time. 

In simulations where the effects of rotation are neglected, we find that material flows along the 
rising tube away from the apex. In the presence of rotation, this type of "draining" flow manifests 
itself as an asymmetry in the axial flow induced by the Coriolis force: the retrograde flow is stronger 
in the trailing side of the loop than it is in the leading side. As the Coriolis force acts upon these 
asymmetric axial flows, and as the buoyant, magnetized plasma attempts to preserve its angular 
momentum, the morphological asymmetry predicted by Moreno-Insertis et al. (1994) and apparent 
in Figure 2 quickly develops. 

Figure 4 shows the transverse components of the Coriolis force, and the flow field in the apex 
cross-section of run MFLL at t = 5.5 (the same time, in naturalized units, as the cross-sections 
of Figure 1). In this run, i?^ = 1, and the domain is again centered at 15° latitude. This implies 
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that Bo = 5.35 x 10^ G, va = 3.37 x 10^ cm s'^, and the unit of time is 2.06 days {t = 5.5 thus 
corresponds to 11.3 days). In an carHer snapshot of run MFLL (t = 1.5), the morphology and flow 
patterns are very similar to those shown for run LFLL in Figure 3. However, as the runs progress, 
differences in secondary and tertiary flows become more pronounced, and the loops begin to evolve 
in a dissimilar manner. This is evident when comparing the first frame of Figure 4 with the frames 
of Figure 1, which highlight the striking differences in morphology that ultimately develop between 
simulations with different values of the magnetic Rossby number. These figures show that as the 
initial axial magnetic field strength of the tube is increased, the tendency for the tube to fragment 
as it rises toward the photosphere is also increased. 



3.1. An Approximate Model 

To better understand the relationship between magnetic Rossby number, latitude of flux emer- 
gence, and loop fragmentation, we approximate the forcing terms of the momentum equation [equa- 
tion (1)] by constant, order-of-magnitude estimates evaluated for our initial conditions, and valid 
for the initial stages of evolution. We then have a simplified expression for the velocity field that 
admits to an analytic solution: 

^ = -2(J2xv) + Q. (5) 

It is important to recognize that this simple approximation will not capture the factor of two 
reduction in acceleration that occurs in the non-rotating case because of "added inertia" effects (see 
Longcope et al. 1996). Nevertheless, it does capture the essence of the dynamic effects introduced 
by the Coriolis force. Let {qxiQy,Qz) represent the Cartesian components of the constant driving 
term Q, and let 9 represent solar latitude. Then the general solution to the above equation is given 

by 

v^{t) = ^sm2nt-vP cos2nt + vP (6) 
Vy (t) = ^ sin 6' cos 2^^^ + sin 9 sin 2^ 

+ {qy - 2nvP sin e)t - ^ sin (7) 
Vz{t) = — ^ COS 9 cos 2Qt — v^^^ cos 9 sin 2^lt 

+ {q^ + 2Q.vP cos 9)t + ^cos9, (8) 

where 

vP = ^{qySm9-q,cos9). (9) 

The total circulation over a cross-sectional slice of the flux tube is given by F = J^u -hs dS = 
§q \ ■ dJt, where C is the closed circuit bounding the cross-sectional surface S, ns is the outward 
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normal to S, and aj=Vxv is the vorticity vector. As the tube begins to rise, oppositely directed 
circulation patterns form on either side of a line defined by the Frenet normal. The total circulation 
of these regions (at the loop apex) is given hy T = <f^,v ■ di ^ dia ~ a < >, if we assume 
that the primary contribution to F stems from v • d£ across the diameter of the tube, a. Since Vz 
is roughly constant across the interior of the tube, we can approximate the circulation (and thus 
the forces acting to pull the tube apart) with our analytic solution for the velocity. Initially, the 
motion of the tube is determined by the buoyancy force, which acts in the z direction. Thus, we 
assume that qx and qy are zero, and express the circulation as: 

^, , avj^ ( .9^ sin(2Qt) , ^ 

We plot our simple approximation of V{t) for four representative values of in Figure 5. 

In order to compare our approximation with the MHD solution, we must first obtain F(t) 
directly from the simulation data. To accomplish this, we follow the methodology described in 
Section 3 of AFF, and first calculate the path of the magnetic flux tube through the background 
plasma. The position of the loop (and its fragments, if applicable) is specified by the magnetic field 
weighted first moments of position along a series of two-dimensional slices that closely correspond 
to the tube's cross-sectional plane. Once the path of the tube is determined (see Section 3 of AFF), 
we construct the Frenet basis vectors, and define two regions on either side of a line defined by the 
Frenet normal where the oppositely directed circulation patterns will develop. From the simulation 
results, we then compute the total circulation across each of these regions. Figure 6 shows V{i) for 
each simulation. One can compare these results with the simple approximation shown in Figure 5. 
It is clear that the general pattern of time dependence of the circulation is successfully reproduced 
by our simple model, though there are significant differences in the details. 

The oscillations in r(t) arc inertia waves (see Batchelor 1967) with a period of t = 2TxR\y (here, 
t is expressed in normalized units of Hy/vA)- Thus, the tendency for the loops to fragment is 
significantly reduced after the circulation begins to fall off at i = ■kRi,/2. The overall evolution 
of the loop depends strongly on the dynamics during the initial portion of its rise, since the early 
differences in secondary and tertiary flows induced by the Coriolis force have a significant impact 
on the eventual evolution of the tube. We find that if the tube is able to rise several of its 
diameters before t reaches 7r-Rf,/2, then the resulting fi-loop will ascend toward the photosphere in 
a fragmented state — otherwise, the tube evolves cohesively. 

We can infer from our approximate model that, on average, the rise speed of the tube slowly 
increases over time. This prediction is confirmed by the simulations, which (for sufficiently high 
choices of and Rm) show the initial slow, linear increase in the rise speed of the tube, modulated 
by inertial oscillations. This behavior persists in the simulations until t becomes comparable to the 
diffusive timescale. If this occurs, magnetic diffusion begins to (artificially) reduce the buoyancy 
of the tube, eventually preventing it from approaching the upper boundary (this occurred in runs 
LFLL and LFHL after t 34.0). Of course, this problem can be mitigated by reducing the 
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coefHcient of magnetic difFusivity in the simulations, but this requires a significant increase in grid 
resolution (to avoid numerical instabilities); and since our analysis concentrates on the early stages 
of loop evolution, such a computationally expensive step is unwarranted. 

3.2. Comparison with Two-dimensional Models 

To investigate the role that three-dimensional geometry plays in the tube's evolution, we 
compare our runs with simulations of buoyant, axially symmetric magnetic flux tubes moving in 
two dimensions. Figure 7 shows the tube cross-section of a 2-D run whose parameters are that 
of the 3-D run MFLL shown in the first frame of Figure 4. This 2-D case can be thought of as 
an infinitely- long, buoyant, axially symmetric fi-loop of zero apex curvature. A comparison of the 
apex cross-section shown in Figure 4 with the cross-section of Figure 7 reveals that a 2-D geometric 
constraint imposed on the solution (even in the early stages of loop evolution) significantly alters the 
dynamic evolution of the tube. This suggests that in the presence of solar rotation, the geometry of 
the J7-loop still plays an important role in determining the morphology of the emerging magnetic 
field — a result consistent with the findings of AFF for similar runs where solar rotation was 
neglected. 

In 2-D, the Coriolis force still acts to suppress the fragmentation of the tube. However, 
the hydrodynamic "lift" force that acts to fragment the tube is overestimated due to the tube's 
effectively infinite axial length scale. It is the binormal component of this force that acts to pull the 
tube apart (see Longcope et al. 1996; Fan et al. 1998; AFF), thus the overall tendency for the tube 
to retain its cohesion during its ascent is reduced as compared to the 3-D cases. In addition, the 
amount that the tube is deflected toward the pole is significantly overestimated in a 2-D geometry. 
This is because in 3-D, the component of the Coriolis force that acts to push the tube poleward 
acts only on the portion of the fi-loop where there arc strong, mainly axial flows (eg. near the 
apex). In general, our calculations support the conclusions of Wissink et al. (2000a) who used 2-D 
simulations to demonstrate that in the presence of the Coriolis force, flux tubes remained somewhat 
more cohesive during their ascent. We note, however, that in 3-D, this effect is more dramatic. 

3.3. Morphology of the Emerging Magnetic Field 

There are several reasons why we must be cautious when making predictions about the be- 
havior of emerging magnetic flux based on data obtained from our simulations. First, the anelastic 
approximation becomes marginal as the flux tube approaches the photospheric boundary, where 

densities are lower, and the acoustic Mach number of the flow approaches unity. Second, we in- 
tended for these simulations to investigate the effects of the Coriolis force alone on the evolution of 
our 0-loop. Thus, we did not include any field line twist along our initial tube, nor did we consider 
the effects of a tube rising through a dynamically convecting background state. Third, our models 
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do not account for radiative losses that occur at the surface. However, given these caveats, we feel 
that we can still make some general predictions regarding expected magnetic field and velocity flow 
patterns if the magnetic flux emerges through a relatively uncluttered region on the solar disk. 

Figure 8 shows the radial component of the magnetic field along horizontal slices taken under- 
neath the apex of the rising i7-loops of runs HFLL, and SSO. This figure emphasizes the differences 
in morphology of emerging magnetic field between a simulation that includes the effects of solar 
rotation (run HFLL), and one that does not (run SSO). In the figures, outwardly directed field is 
displayed in white, and inwardly directed field is displayed in black. Superimposed over the each 
grey-scale image is the transverse velocity field. One of the most striking effects of the Coriolis force 
is the strong flow that rapidly develops in both the magnetized and unmagnetized plasma directed 
from the leading toward the trailing polarity (the leading polarity moves in the direction of solar 
rotation as the loop continues to rise). This flow is simply the result of the Coriolis force acting on 
the rising fluid just below the apex of each loop. This "LTT" (Leading To Trailing polarity) flow 
is evident in all of our simulations except in cases where solar rotation is neglected. 

The Coriolis force manifests itself in several other interesting ways. First, since the trailing 
edge of the loop emerges at a steeper angle with respect to the surface than the leading edge (an 
effect most pronounced in the runs where Rb is small), the leading polarity becomes "spread-out" 
in the east-west direction (along x). Second, in the cases where is small, we observe that the 
leading polarity is positioned slightly closer to the equator than the trailing polarity — behavior 
that is roughly consistent with the "Joy's Law" result from thin flux tube models (see D'Silva & 
Choudhuri 1993; Fan et al. 1994; Schiissler et al. 1994; Cahgari et al. 1995; Fan & Fisher 1996; 
Fisher et al. 2000). Third, the north-south dispersion of the emerging flux is greatly reduced when 
solar rotation is included in the simulations. This reflects the increased cohesion of the flux rope 
during its ascent through our model convection zone. 

4. Summary and Conclusions 

AFF found that the degree of fragmentation of a rising O-loop depends upon the 3-D geometry 
of the loop — the greater the apex curvature, the lesser the degree of fragmentation for a fixed 
amount of initial field-line twist. We are now able to extend this analysis to a rotating model 
convection zone, and investigate how the degree of apex fragmentation depends on the axial field 
strength of the initial magnetic flux tube. Assuming a given rotation rate, we find that: 

1. A magnetic flux tube that rises buoyantly toward the photosphere in a rotating convection 
zone rises more slowly, and is better able to retain its cohesion than an identical flux tube 
that rises through a non-rotating convection zone — even in the absence of initial field line 
twist. If the tube is sufficiently buoyant so that it is able to rise a distance of several tube 
diameters before t = 7r/4J7 (where O is an estimate of the relevant solar rotation rate at a 
given latitude), then the tube will fragment during its rise toward the photosphere. 
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2. In the absence of forces due to convective motions, the stronger the initial axial field strength 
of the tube, the more fragmentation occurs as the tube rises. 

3. An approximate model is able to predict the initial evolution of the circulation, and hence can 
effectively characterize the formation and subsequent interaction of the oppositely-directed 
vortex pairs of a fragmented S7-loop. 

4. One must be careful not to over-interpret results from two-dimensional simulations of buoyant 
magnetic flux tubes. In a two-dimensional geometry, hydrodynamic forces due to vortex 
interaction are overestimated. The amount of poleward deflection due to the Coriolis force is 
also overestimated. 

5. If magnetic flux emerges as an fi-loop, the Coriolis force will induce a relatively strong flow 
that is directed from the region of leading polarity to the region of trailing polarity. This flow 
is present in both the magnetized and unmagnetized plasma. 
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Fig. 1. — Apex cross-sections of |B| for a simulation with Rf, = 0.5 at 15° latitude (left frame, run 
LFLL), and a simulation neglecting solar rotation (right frame; labeled run SSO to be consistent 
with the naming convention of AFF). |B| is given in terms of the field strength along the axis of 
the initial tube, and t is given in units of the ratio of the pressure scale height at the base of the 
domain to the Alfven speed along the axis of the initial tube. The coordinate axes are given in 
terms of computational zones, though only a portion of the domain is shown. 
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Fig. 2. — A volume rendering of |B| for the Jl-loop of run HFLL at t 
computational domain is shown. 



= 8.53. The entire rectangular 
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LFLL APEX t=1 .5 




Fig. 3. — Apex cross-sections of |B| along with the transverse Coriohs force (left frame), and the 
corresponding velocity field (right frame) for run LFLL (Rf, = 0.5 at 15° latitude) at t = 1.5. 
Quantities are given in naturalized units (see text for details). In the right frame, axial flows 
depicted in red indicate velocities receding from the viewer, while blue colors indicate flow moving 
toward the viewer. Here, max|va;|/max|vT| ~ 1-5; that is, the axial flows are stronger than the 
transverse flows. The angular velocity vector is displayed in the lower left-hand corner of the left 
frame, and applies to the entire figure. 



- 20 - 




Fig. 4. — Same as Figure 3 for run MFLL (Ri^ = 1 at 15° latitude) at t = 5.5 (the same time as 
displayed in Figure 1). Here, nicixl I /majxl I ^ 1.1. 
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Fig. 5. — Analytic solution for the total circulation of a tube "fragment" as a function of time for 
four representative values of the magnetic Rossby number. The quantities are given in naturalized 
units. 



- 22 - 



APEX (simulations) 




LFHL 
LFLL 



1 2 3 4 5 6 



Fig. 6. — 3-D MHD solution for the total circulation of a tube "fragment" at latitudes of 15° and 
30° as a function of time and magnetic Rossby number. The quantities are given in naturalized 
units. 
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Fig. 7. — A cross-section of |B| at t = 5.5 (in normalized units) for a 2-D run that begins with the 
initial magnetic flux tube of run MFLL. 
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Fig. 8. — A grey-scale image of the radial magnetic field for a horizontal slice just below the apex 
of the fi-loops of the final timesteps of runs HFLL (top frame) and SSO (bottom frame). The 
vector field represents the transverse (x-y) flow field. The x and y scales are given in terms of 
computational zones. The grid is uniform; each zone is roughly 1 Mm in size. There is a very slight 
tilt of approximately 2.5 degrees from horizontal between the center of the bipolar regions. 
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Table 1. Physical Parameters 



Label 


Resolution^ 




Bo{G) 


Rf, 




H,!^ (cgs) 


LFHL 


256x128x128 


30 (13.5135) 


2.5966x10'' 


0.5 


1700.07 


5.7805x10^0 


LFLL 


256x128x128 


15 (13.9104) 


2.6728x10^ 


0.5 


1750.00 


5.7805x10^° 


MFHL 


256x128x128 


30 (13.5135) 


5.1932x10^ 


1.0 


3400.15 


5.7805x10^° 


MFLL 


256x128x128 


15 (13.9104) 


5.3457x10^ 


1.0 


3500.00 


5.7805x10^° 


HFHL 


256x128x128 


30 (13.5135) 


1.0386x10^ 


2.0 


3400.15 


1.1561x10^^ 


HFLL 


256x128x128 


15 (13.9104) 


1.0691x10^ 


2.0 


3500.00 


1.1561x10^^ 



''Latitude, degrees 

^Observed rotation rate, degrees/day — data from MDI 2dRLS inversion (Scherrer 
1997) 



